Shock-Like Dynamics of Inelastic Gases 
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We provide a simple physical picture which suggests that the asymptotic dynamics of inelastic 
gases in one dimension is independent of the degree of inelasticity. Statistical characteristics, includ- 
ing velocity fluctuations and the velocity distribution are identical to those of a perfectly inelastic 
sticky gas, which in turn is described by the inviscid Burgers equation. Asymptotic predictions of 
this continuum theory, including the t~ 2 ^ 3 temperature decay and the development of discontinuities 
in the velocity profile, are verified numerically for inelastic gases. 
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Gases of inelastically colliding_particles model the dy- 
namics of granular materials geophysical flows [3|, 
and large-scale structure of matter in the universe 
Typically, a fraction of the kinetic energy is dissipated 
in each collision, leading to interparticle velocity correla- 
tions, a clustering instability and in the absence of 
external energy input, an inelastic collapse JlO|-p^|. The 
last feature presents an obstacle to long-time simulations, 
as an infinite number of collisions occur in a finite time. 

In this Letter, we propose that a freely evolving in- 
elastic gas is asymptotically in the universality class of 
a completely inelastic, sticky gas. Specifically, the tem- 
perature decreases in time as i~ 2 over an intermediate 
range, but asymptotically decays as i~ 2 / 3 . To test this 
hypothesis, we employ a simulation in which collisions 
between particles with sufficiently small relative veloci- 
ties are perfectly elastic. This method allows us to bypass 
the inelastic collapse and probe the asymptotic regime. 

We consider N identical point particles undergoing in- 
elastic collisions in a one-dimensional periodic system of 
length L. The particles have typical interparticle spac- 
ing xq = L/N and their typical velocity is vq. We em- 
ploy dimcnsionless space and time variables, x — > x/xq, 
and t — ► tvo/xo, thereby rescaling the ring length to N. 
Inelastic and momentum conserving collisions are imple- 
mented by changing the sign of the relative velocity and 
reducing its magnitude by a factor r — 1 — 2e, with 
< r < 1, after each collision. It is convenient to view 
the particle identities as "exchanged" upon collision, so 
that in a perfectly elastic collision the particles merely 
pass through each other, while for a small inelasticity 
each particle suffers a small deflection. The outcome of a 
collision between a particle with velocity v and another 
particle with velocity u is therefore 



v — e(v — u). 



(1) 



The granular temperature, or velocity fluctuation, 
T(t) = (v 2 (t)} — (u(t)) 2 , can be estimated in the in- 
termediate time regime by considering the outcome of 
a single collision under the assumption that the system 
remains homogeneous. In each such collision, the en- 
ergy lost is AT oc —e(Av) 2 , with Av the relative veloc- 
ity, while the time between collisions is lj Av. Assum- 
ing homogeneity, we neglect fluctuations in the mean- 



free path £ = 1 and posit a single velocity scale so that 
v ~ Av ~ T 1 / 2 . The temperature therefore obeys the 
rate equation dT/dt cx — eT 3 / 2 



giving 



T{t) ~ (1 + Aety 



(2) 



with A a constant of order unity M. For small times 



e dissipation is negligible and the tem- 



perature does not evolve - the gas is effectively elastic. 
For larger times, the dissipation leads to a e~ 2 t~ 2 tem- 
perature decay. 

However, this behavior cannot be valid asymptoti- 
cally, as the temperature must decrease monotonically 
with increasing dissipation. Moreover, the temperature 
is bounded from below by that of the perfectly inelas- 
tic gas with a vanishing restitution coefficient, r = 0. 
For such a sticky gas, the temperature decays as i~ 2 / 3 
and the typical cluster mass grows as t 2 / 3 @. This be- 
havior is reminiscent of diffusion-controlled two species 
annihilation, where a small reaction probability results 
in a homogeneous intermediate time regime where the 
density follows aT 1 mean-field decay, even for low spa- 
tial dimension d. However, at long times single-species 
domains which are opaque to opposite-species particles 
form and a slower t -d / 4 density decay follows [ p~7| . 

For the inelastic gas, we argue that the role of the re- 
action probability is played by e. For small e, a particle 
can penetrate through a domain of N < N c (e) ~ e _1 
coherently-moving particles without experiencing a sub- 
stantial deflection. The critical cluster size N c (e) may be 
estimated by considering a collision between a moving 
particle and a cluster of N stationary particles. From 
Eq. (0), each collision between the incident particle and 
the next particle in the cluster reduces the incident parti- 
cle velocity by roughly e. After A" collisions the incident 
velocity is v^ ss 1 — Ne. For the particle to pass through 
the cluster, the number of particles must therefore be less 
than e _1 . It is in this range of cluster sizes that the sys- 
tem remains spatially homogeneous and the mean-field 
decay T ~ e^ 2 i~ 2 holds. 

However, once the cluster size is larger than A r c (e) ~ 
e _1 , an incident particle is "absorbed" and the decay 
follows that of the perfectly inelastic gas. That is, do- 
mains larger than N c (e) are opaque and present an effec- 
tive restitution coefficient r e g = to incident particles. 
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We argue that a similar sticking mechanism also gov- 
erns cluster-cluster collisions. The crossover time i s tick 
between these two regimes is obtained by matching the 
intermediate and long-time temperature decays, e~ 2 t~ 2 
and t~ 2 / 3 , to give t stick ~ e~ 3/2 . 

These arguments suggest the temperature decay 

!1 t <C = idissip! 

e- 2 t- 2 e- 1 < t < e~ 3/2 = t stick ; /o\ 
£-2/3 e -3/a « f JV3/2. W 
^3/2 < f 

The last regime reflects the final state of a finite N- 
particle system, namely a single cluster of mass m = N , 
velocity u ~ TV -1 / 2 , and therefore energy T ~ v 2 ~ iV -1 . 
This final velocity follows from momentum conservation 
in which the total momentum P is the sum of N indi- 
vidual random momenta of order unity. Consequently, 
P oc N 1 / 2 and v = P/m ~ N' 1 / 2 . 

The above crossover picture applies equally well to 
moderately inelastic gases where both idissip and i s ti c k 
are of order unity and the asymptotic behavior is real- 
ized immediately. While weakly inelastic systems with a 
small number of particles N < e _1 will avoid the cluster- 
ing regime and follow the t~ 2 cooling law indefinitely, the 
t~ 2/3 sticky gas regime is always reached in the thermo- 
dynamics limit, N — > co. Therefore, the t~ d / 2 decay con- 
jectured in |l8|,|l9| based on two- and three-dimensional 
simulations does not extend to lower dimensions. 

To probe the long-time behavior, we performed numer- 
ical simulations of N particles which are initially equally 
spaced (Ax = 1) and uniformly distributed (on [—1,1]) 
in velocity. We implemented an event-driven simula- 
tion, keeping the collision times always sorted to facilitate 
identification of the next event. To circumvent the inelas- 
tic collapse, elastic collisions were implemented whenever 
the relative velocity of the colliding particles fell below a 
pre-specified threshold, Av < S ||J2(|- In fact, the resti- 
tution coefficient for deformable spheres does approach 
unity when Av — > as a consequence of the nonlinear 
Hertz contact law f2l}| . 

Fig. 1 shows that the temperature of the freely cool- 
ing inelastic gas asymptotically decays as t~ 2 / 3 , indepen- 
dent of the restitution coefficient. Moreover, the time 
scale over which this decay occurs diverges in the limit of 
vanishing dissipation. We also simulated the completely 
inelastic gas (r = 0) where particles aggregate (and con- 
serve momentum) upon collision. This gives an identical 
asymptotic temperature decay as the partially inelastic 
gas. From Fig. 1, notice that the two crossover times 
^dissip ~ e_1 and istick ~ e -3 / 2 are consistent with the 
data for the cases r = 0.9 and 0.99, and that the tem- 
perature is of the appropriate order T(£ st i c k) ^ e at the 
homogeneous-sticky crossover. On the other hand for 
r = 0.5, the intermediate t~ 2 regime no longer exists and 
only the sticky gas behavior is realized. 
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FIG. 1. Temperature T(t) versus time for the freely cool- 
ing inelastic gas with restitution coefficients r = 0.5, 0.9, and 
0.99. The data represents averages over 10 realizations of 
N = 10 6 particles and 10 4 collisions per particle. Also shown 
is a simulation for a sticky gas (r = 0). A dashed line of 
slope —2/3 is plotted as reference. Least-square fits to the 
post-crossover data with velocity threshold 8 = 10~ 2 yield 
the decay exponents 0.67, 0.67, and 0.66 for r — 0, 0.5, and 
0.9 respectively. 
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FIG. 2. Role of threshold velocity and the collision mecha- 
nism. Temperature decay for r = 0.9, N = 10 6 with different 
threshold velocities S — 10 -2 , 8 — 10~ 3 , and different sub- 
threshold collision mechanisms (both sticky and elastic). 

The crossover between the homogeneous gas regime 
and the sticky gas regime is ultimately related to the in- 
elastic collapse. In Fig. 1, simulation results for purely 
inelastic collisions are also shown (dots); these coincide 
with the results for velocity threshold 6 = 10~ 2 . How- 
ever, as soon as the first cluster is formed, the number 
of collisions becomes infinite and the simulation does not 
proceed any further in time. Indeed, the last data point 
for these simulations marks the transition to the sticky 
gas regime. This provides an additional confirmation for 
the crossover picture, as the typical cluster mass of the 

sticky gas, m(4tick) ~ ^stick ~ e_1 ' matches the critical 
mass for the appearance of opaque clusters. 



2 




X 

FIG. 3. Space-time evolution of a 500 particle system with 
r = 0.9 and 8 = 10~ 2 , up to t = 600. 

To validate the simulation method, we checked that 
results are independent of the cutoff value (provided it 
is sufficiently small) as well as the sub-threshold colli- 
sion mechanism (Fig. 2). In principle, the results can be 
trusted as long as the typical velocity is much larger than 
the cutoff, v ~ i -1 / 3 3> S, i.e., up to time i va iid ~ S~ 3 - 
As shown in the figure, the results for 5 = 1CP 2 and 
1CP 3 nearly coincide until t = 10 6 , consistent with our 
expectation. Furthermore, the space-time evolution of a 
weakly inelastic gas illustrates how aggregation eventu- 
ally dominates (Fig. 3). 
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FIG. 4. Scaling of the velocity distribution for restitution 
coefficients r = 0, 0.5, and 0.9, at times t = 10 3 , 10 4 , 10 5 . 
The nine data sets represent averages over 200 realizations in 
a system of N — 10 5 particles and cutoff S = 10~ 3 . 

We now investigate whether the velocity distribution, 
and not merely the overall velocity scale, is also indepen- 



dent of r. We therefore computed this distribution for 
r = 0.9 (weakly inelastic), r — 0.5 (moderately inelas- 
tic), and r = (perfectly inelastic) at three very differ- 
ent times which are well into the clustering regime. For 
the r = case, the cluster velocity was weighted by the 
cluster mass, to compare with the r > cases. As shown 
in Fig. 4, the normalized velocity distribution 

P(v,t) ~ t 1 / 3 ^ 1 / 3 ), (4) 

is described by an identical scaling function $(z) for these 
widely different values of r. This universality provides 
further confirmation that the asymptotic behavior for 
any r < 1 is governed by the r = "fixed point" . 

Further insights about the behavior of the inelastic gas 
are provided by the connection to the Burgers equation 

P, p2[ . Since sticky gases are described by the inviscid 
— » 0) limit of the Burgers equation 

Vt + vv x = vv xx , (5) 

supplemented by the continuity equation p t + {pv) x = 0, 
we conclude that this continuum theory also describes 
the asymptotics of the inelastic gas in the thermody- 
namic limit. The Burgers equation may be reduced to 
the diffusion equation by the Hopf-Cole transformation 
v = — 2f(lnu) x , and therefore is solvable. In our case, 
the relevant initial condition is delta-correlated velocities 
(vq(x)vq(x')) = S(x — x'). The resulting velocity profile 
is discontinuous, and the corresponding shocks can be 
identified with clusters in the sticky gas. Indeed, both 
shock coalescence processes and cluster-cluster collisions 
in the sticky gas conserve mass and momentum. 

The relation to the Burgers equation is useful in sev- 
eral ways. First, statistical properties of the shock coa- 
lescence process have been established analytically [ p3| . 
For example, the tail of the particle velocity distribution 
([|) is suppressed according to 

$(z) ~ exp(-const. x \z\ 3 ), \z\ > 1. (6) 

This behavior can be understood by considering the den- 
sity of the fastest (order unit velocity) particles. For 
such a particle to maintain its velocity to time t, it must 
avoid collisions. This requires that an interval of length 
oc t ahead of the particle must be initially empty p4[ . 
For an initially random spatial distribution, the prob- 
ability of finding such an interval decays exponentially 
with length; thus P(l,t) ~ exp(— const, x t). Using 
$(z) ~ exp(— const, x |z| 7 ) and z — vt 1 / 3 then yields 
7 = 3. Interestingly, over most of the range of scaled 
velocities, the numerically obtained velocity distribution 
deviates only slightly from a Gaussian, reflecting the 
small constant in (^|) p3[ . 

Another important prediction of Eq. (|^) is that the 
velocity is linear in the Eulerian coordinate x and the 
Lagrangian coordinate q(x,t) 

I + \ x-q{x,t) 
v(x,t) = . (7) 
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This form also characterizes the asymptotic velocity pro- 
file of inelastic gases. Fig. 5 shows such a sawtooth ve- 
locity profile from an inelastic gas simulation. The slopes 
of the linear segments of the profile are consistent with 
the i _1 prediction of Eq. (Q). The inelastic collapse is 
simply a finite time singularity characterized by the de- 
velopment of a discontinuity in the velocity profile, i.e., 
a shock. 
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FIG. 5. Shock profile of inelastic gases. Density and veloc- 
ity are plotted at time t = 10° in a system with N = 2 x 10 
particles, r = 0.99 and 8 — 10 -4 . A line of slope t^ 1 is plot- 
ted for reference. The number of shocks is consistent with the 
expected number Nt~ 2 ' 3 = 9. 

In higher dimensions as well, the temperature of an 
inelastic gas is a monotonically increasing function of 
r and hence, it is bounded from below by the r = 
case. Therefore, we speculate that r = remains the 
fixed point in higher dimensions. On the other hand, 
the Burgers equation vt + v • Vv = i/V 2 v approximately 
describes the sticky gas in the limit v — > The 
known t~ d > 2 temperature decay of the Burgers equation 
j| , valid for 2 < d < 4 (with possible logarithmic correc- 
tions at the crossover dimensions), then yields 



T(t) 



1 i«( 

-2+-2 r -l 



e H 

t -d/2 

N- 1 



e ^ <t <e- 4 /( 4 - d ); 
e -4/(4-d) <<; t <<; N 2/d. 

N 2 ' d < t. 



Interestingly, both the decay exponents |ll|[l9|, the for- 
mation of string- like clusters |5 Hll] , |20|| , and even the pos- 
sibility of a percolating network of clusters p5| , fea- 
tures that were found primarily numerically, are all pre- 
dicted by the Burgers equation. Additionally, the criti- 
cal cluster size increases with the dimension according to 
, suggesting that the inelastic collapse 
= 4, and that the homogeneous 
2 holds indefinitely above this 



Nr. 



-2d/(4-d) 



is avoided when d > d c 



gas behavior T ~ e~H~ 
critical dimension. 

In summary, our results suggest that the asymptotic 



behavior of a one-dimensional inelastic gas with many 
particles is governed by the r = sticky gas fixed point, 
and that the appropriate continuum theory is the invis- 
cid Burgers equation. This connection provides several 
exact statistical properties of inelastic gases. Conversely, 
inelastic gases may provide a useful tool to study shock 
dynamics. The suggestive behavior of the inelastic gas in 
high dimensions deserves careful investigation. 
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